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The circular and linear magnetic birefringences corresponding to the Faraday and the Cotton-Mouton effects, 
respectively, have been measured in xenon at A = 1064 nm. The experimental setup is based on time dependent 
magnetic fields and a high finesse Fabry-Perot cavity. Our value of the Faraday effect is the first measurement 
at this wavelength. It is compared to theoretical predictions. Our uncertainty of a few percent yields an 
agreement at better than Icr with the computational estimate when relativistic effects are taken into account. 
Concerning the Cotton-Mouton effect, our measurement, the second ever published at A = 1064 nm, agrees 
at better than Icr with theoretical predictions. We also compare our error budget with those established for 
other experimental published values. 
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I. INTRODUCTION 

Magnetic birefringence corresponds to an anisotropy 
of the (generally complex) refractive index induced in a 
medium by a magnetic field.— A circular birefringence 
arises when the magnetic field changes the angular ve¬ 
locity of the two eigen modes of polarization in which a 
linearly polarized beam is split, without deforming them. 
The net result is a rotation of the plane of linear polar¬ 
ization, a phenomenon seen also in absence of external 
fields in chiral samples (natural optical rotation). When 
the presence of the external magnetic field yields a differ¬ 
ent phase of two perpendicular components of the linear 
polarization vector, the net result is the appearance of 
an ellipticity, and we are observing an example of linear 
birefringence. 

Two well known examples of magnetic birefringences 
are the Faraday and the Cotton-Mouton effects. The 
former corresponds to a circular birefringence induced by 
a longitudinal magnetic field il|| (aligned parallel to the 
direction of propagation of light). After going through 
the birefringent medium, the real part of the index of 
refraction for left circularly polarized light n_ is different 
from that for right circularly polarized light n+. The 
difference Anp = n_ — n+ is proportional to i?|| 

Anp = ( 1 ) 

fcp being the circular magnetic birefringence per unit 
magnetic field intensity. For historical reason, the Fara- 
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day effect is usually given in terms of the Verdet constant 


where A is the light wavelength. On the other hand, the 
Cotton-Mouton effect corresponds to a linear magnetic 
birefringence induced by a transverse magnetic field B ±. 
The field induces a difference between the real parts of 
the refraction index for light polarized parallel with re¬ 
spect to that polarized perpendicular to the magnetic 
field. The difference AncM = ^|| ~ ^-t is proportional to 
the square of the magnetic field 

AncM = (3) 

with fccM the linear magnetic birefringence per square 
unit magnetic field intensity. 

For the Cotton-Mouton effect, fccM has two contribu¬ 
tions, the first one due to the distortion of the electronic 
structure while the second one corresponds to a partial 
orientation of the molecules. When working in the condi¬ 
tions of constant volume, the orientational contribution 
is proportional to the inverse of the temperature T, and 
it usually dominates, often hiding the first temperature 
independent contribution. For axial molecules, for exam¬ 
ples, fccM is given by the expression^ 

Above Na is the Avogadro constant, Vm the molar vol¬ 
ume, fee the Boltzmann constant, cq the electric con¬ 
stant, A 77 the frequency dependent hypermagnetizabil- 
ity anisotropy, Aa the optical electric dipole polariz¬ 
ability anisotropy, and Ay the magnetic susceptibility 
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anisotropy. For spherical molecules or for atoms, such 
as xenon, however, the temperature dependent contribu¬ 
tion vanishes. Measurements on noble gases, for example, 
allow to focus on the hypermagnetizability anisotropy 
Ar] term. On the other hand, since the Langevin-type 
orientational term vanishes, the magnetic birefringence 
is much lower than the one observed in non spherical 
molecules. From an experimental point of view, measure¬ 
ments on such gases require a very sensitive apparatus, 
with a AncM of the order of 10 “^® for helium and 10 “^"^ 
for xenon at one atmosphere and with a magnetic field 
of one Tesla. In comparison, Anp is typically 10® bigger. 


The computational determination of the Verdet con¬ 
stant and of the Cotton-Mouton effect requires the far- 
from-trivial calculation of higher-order response func¬ 
tions,—!^ and it has often served as test bed for the vali¬ 
dation of new electronic structure methods. For atoms, 
in order to obtain accurate results one must properly 
account for the appropriate description of one-electron 
(basis set), N-electron (correlation) and relativistic ef¬ 
fects. As far as correlation is concerned, coupled cluster 
(CC) methods are nowadays among the most accurate 
tools in electronic structure theory— Both birefringences 
treated here, and in particular the Cotton-Mouton effect, 
require a good description of the outer valence space of 
the system at hand, and therefore the presence of diffuse 
functions in the one-electron basis set is mandatoryi^ 
Whereas for light atoms relativistic corrections are mi¬ 
nor, their importance increases and they become signif¬ 
icant for heavier atoms. For example, Ekstrom et al^ 
have calculated that for helium the relativistic effects add 
—0.03% to the non-relativistic Verdet value. For xenon, 
the heaviest non-radioactive noble atom, relativistic cor¬ 
rections add 3 to 4%, depending on the chosen wave¬ 
length. In this case, relativistic effects cannot be ignored 
in accurate calculations. 


In this article, we report both measurements and cal¬ 
culations of Faraday and Cotton-Mouton effects at A = 
1064 nm. We perform the first measurement of the Fara¬ 
day effect of xenon at this wavelength, and our esti¬ 
mate bears an uncertainty of a few percent. Concerning 
the Cotton-Mouton effect, our measurement, the second 
ever published at A = 1064 nm, agrees at better than Ict 
with theoretical predictions and we also compare our er¬ 
ror budget with those established for other experimental 
published values. Our theoretical predictions, that can 
be considered of state-of-the-art quality, were obtained at 
the coupled cluster singles and doubles (CCSD)!^— and 
coupled cluster singles, doubles and approximate triples 
(CC3)^“— levels of theory, and they include estimates 
of relativistic effects. For both effects, our theoretical 
predictions are within Icr of our experimental data. 


II. EXPERIMENTAL SETUP 

A. Principle of the measurement 

Experimentally, we determine the Faraday and the 
Cotton-Mouton effects by measuring, respectively, the ro¬ 
tation induced by a longitudinal magnetic field and the 
ellipticity induced by a transverse magnetic field on an 
incident linear polarization. For small angles, the in¬ 
duced rotation 9^ depends on the circular birefringence 
as follows 

6 >f = tt^Auf, (5) 

where Lb is the length of the magnetic field region. The 
induced ellipticity ipcM is related to the linear birefrin¬ 
gence by the formula: 

Lb 

ipCM = TT—AncMsin20p, (6) 

where 0p is the angle between the light polarization and 
the magnetic field. 

B. General setup 

The apparatus has already been described in detail 
elsewhere i^^d^ Briefly, light comes from a Nd:YAG laser 
at A = 1064nm (see Fig.[T]). It is linearly polarized by 
a first polarizer P, before going through a transverse or 
a longitudinal magnetic field. The polarization is then 
analyzed by a second polarizer A, crossed at maximum 
extinction compared to P. The beam polarized parallel to 
the incident beam, reflected by the polarizer A as the or¬ 
dinary ray, is collected by the photodiode Pht. Its power 
is denoted by It- The beam polarized perpendicular to 
the incident beam (power Ig), corresponding to the ex¬ 
traordinary ray that passes through the polarizer A, is 
collected by the low noise and high gain photodiode Phg. 


Pht 



FIG. 1. Experimental setup. EOM = electro-optic modula¬ 
tor; AOM = acousto-optic modulator; PDH = Pouird-Drever- 
Hall; Ph = photodiode; P = polarizer; A = analyzer. See text 
for more details. 

This setup has been designed to measure the linear 
magnetic birefringence of vacuun>i^ and its sensitivity 
allows to perform precise measurements on gasest ^®d^ 
All the optical components from A to P are placed in 
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an ultrahigh-vacuum chamber. To perform birefringence 
measurement on gases, we fill the vacuum chamber with 
a high-purity gas. For this particular measurement, we 
have used a bottle of xenon with a global purity higher 
than 99.998%. 


C. Fabry-Perot cavity 

Magnetic birefringence measurements on dilute gases 
are difficult, especially at low pressure, because one has 
to detect very small variations of light polarization. To 
increase the measured signal, one needs high magnetic 
fields. One also needs an as large as possible path length 
in the field Lb (cf. Eqs (l5|) and ([51)). To this end, optical 
cavities are used to trap light in the magnetic field region 
and therefore enhance the signal to be measured. 

As shown in Fig. [I] the cavity is formed by two mir¬ 
rors Ml and M 2 , placed at both sides of the mag¬ 
netic field region. The laser frequency is locked to 
the cavity resonance frequency, using the Pound-Drever- 
Hall techniqueji^ The electro-optic modulator generates 
10 MHz sidebands and the signal reflected by the cavity 
is detected by the photodiode Ph,.. The laser frequency 
is adjusted with the acousto-optic modulator, the piezo¬ 
electric and the Peltier elements of the laser. 

This cavity increases the distance traveled by light in 
the magnetic field by a factor 2F/'k, where F is the cavity 
finesse. Therefore, the rotation induced by the longitu¬ 
dinal magnetic field becomes 

2F 

0 F(t) = —0F(t), (7) 

TT 

with 9f the rotation acquired without any cavity. In 
the same way, the ellipticity induced by the transverse 
magnetic field becomes 

2F 

'I'CM(i) = — ’fpcuit), ( 8 ) 

TT 

with ipcM denoting the ellipticity acquired without any 
cavity. The cavity finesse is inferred from the measure¬ 
ment of the photon lifetime r inside the cavityi^ 

F = 27rA^®’^T, (9) 

with the cavity free spectral range. For the Faraday 
effect, the cavity finesse was about F = 475 000. For the 
Cotton-Mouton effect, two sets of mirrors were used with 
a respective finesse of about 400 000 and 480 000. 


D. Raw signals 

We measure the circular and the linear magnetic bire¬ 
fringence by measuring the ratio L /It 

r W = -I- [F -I- TcmW]^ + [s + 0F(i)]^- (10) 

h,iW 


As said previously, le (A) corresponds to the power of 
light polarized perpendicular (parallel) to the incident 
beam. The subscript f indicates that we need to take 
into account the cavity filtering, as explained in details 
in previous papersi ^^d^ The term corresponds to the 
extinction ratio of polarizers P and A, F is the total static 
ellipticity due to the cavity mirrors and e is the static an¬ 
gle between the major axis of the elliptical polarization 
and the incident polarization. The extinction ratio and 
the static birefringence are measured before each mag¬ 
netic pulse. The static angle e can be estimated but its 
value is not needed for the analysis. 


III. CIRCULAR MAGNETIC BIREFRINGENCE 
A. Magnetic field 

The magnetic field is generated by a solenoid previ¬ 
ously used for Faraday effect measurement in heliumJ^ 
Its characteristics have already been explained in de- 
tailsfi^ Here we just briefly recall its main features. It 
generates a longitudinal magnetic field with an equiva¬ 
lent length Lb = (0.308 ± 0.006) m at Icr. This mag¬ 
netic field is modulated at the frequency iz = 18 Hz: 
H|| = H|| 0 sin(27rjzt-|-^). The rotation of the polarization 
due to the Faraday effect is thus given by 

0F = 00 sin(27r^f -I- ^), (11) 

2F 

with 00 = — yB\\ qBb- (12) 

TT 


B. Data analysis 


Expanding Eq. (flUl) . the raw signal becomes 

.I^=a^+T^ + e^ + 2eQFit) + QUt)- ( 13 ) 

A,f(r) 


This gives three main frequency components: a DC sig¬ 
nal, a signal at the frequency v, and a signal at the double 
frequency 2iy. To measure the Verdet constant, we use 
the amplitude of the signal at 2i*^ 


A 2 ,, — 


02 


2 Vl+(f)' 


(14) 


where Vc = l/47rT is the cavity cutoff frequency, intro¬ 
duced to take into account the cavity filteringJ^ A 21 / is 
measured for different magnetic field amplitudes, from 0 
to about 50 X 10“^ T. The whole is fitted by p. The 

Verdet constant finally depends on the measured experi¬ 
mental parameters as follows 


i/(r,p) 


IKv 

1 -I- (87rrzz)^ 

1/4 

^ 2 

2tAFSRLb 


(15) 


where T and P are respectively the temperature and 
pressure of the gas. 
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TABLE I. Parameters and their respective relative A- and 
B-type uncertainties at Icr that have to be measured to infer 
the value of the Verdet constant V. Typical values are given 
at P = 5 X 10“® atm. 


Parameter 

Typical 

value 

Relative 

A-type 

uncertainty 

Relative 

B-type 

uncertainty 

r [ms] 

1.14 

2.0 X 10-^ 


Kv [rad T-^] 

1.07 

3 X 10-3 

3.2 X 10-^ 

[mHz[ 

65.996 


3 X 10-'‘ 

Lb [m] 

0.308 


1.9 X 10-^ 


V X 10® 1.66 

1.8 X 10-^ 

2.5 X 10-^ 

[radT-3m-3] 




C. Measurement and error budget 

The A- and B-type uncertainties associated to the mea¬ 
surement of V are detailed in Tab. They are given 

at ItT (coverage factor k = 1). The A-type uncertainty is 
dominated by the photon lifetime uncertainty. The main 
contributions of the B-type uncertainty comes from the 
uncertainty of the magnetic length and of the fit constant 
Ky which includes the B-type uncertainty of the mag¬ 
netic field and of the photodiodes conversion factor—. 

We have measured the Verdet constant in xenon at 
T = (294 ± 1)K and for 5 pressures from 1.01 x 10“^ 
to 5.01 X 10“^ atm. In this range of pressure, xenon can 
be considered as an ideal gas and the Verdet constant is 
thus proportional to the pressure. Data are fitted by a 
linear equation: 


I/(T,P) = V"P, (16) 

giving a normalized Verdet constant (P = 1 atm) at A = 
1064 nm and T = (294 ± 1) K 

V*" = (3.31 ± 0.09) X 10-3 atm-^rad T-^m-\ (17) 

The uncertainty is given at Icr and is detailed in Tab.llll 
With a scale law on the gas density, this corresponds to 
a normalized Verdet constant at T = 273.15 K of 

= (3.56 ± 0.10) X 10-3 atm-^rad T-^m-^ (18) 

Using Eq. ([2]), we can also give the normalized Faraday 
constant at T = 273.15 K 

k^ = (1.21 ± 0.03) X 10-^ atm-3T-3. (19) 


IV. LINEAR MAGNETIC BIREFRINGENCE 


TABLE 11. Parameters and their respective relative A- and 
B-type uncertainties at Icr that have to be measured to infer 
the value of the normalized Verdet constant U". The uncer¬ 
tainty given by the linear fit takes into account the A-type 
uncertainty of V. 


Parameter 

Typical 

value 

Relative 

A-type 

uncertainty 

Relative 

B-type 

uncertainty 

V X 10® 
[radT-^m-i] 

1.66 

1.8 X 10-^ 

2.5 X 10-2 

P X 10® 

[atm[ 

5 


2 X 10-® 

linear fit 
xlO® 

[atm-^rad 

T-3m-3] 

3.31 

1.5 X 10-2 


U" X 10® 

3.31 

1.5 X 10-2 

2.5 X 10-2 


[atm ^rad 


National Laboratory (LNCMI-Toulouse, France) for the 
measurement of the vacuum magnetic birefringence. This 
coil has been presented and discussed in great details in 
several previous papers— Very briefly, the magnet de¬ 
livers a pulsed magnetic field over an equivalent length 
Lb of 0.137m. The total duration of the pulse is about 
10 ms with a maximum reached within 2 ms. For the 
present measurements, a maximum magnetic field of 3 T 
has been used. Finally, the high-voltage connections can 
be remotely switched to reverse the direction of the field. 
Thus we can set Bx parallel or antiparallel to the x di¬ 
rection, as shown in Fig.[TJ 


B. Data analysis 


The data analysis follows the one described for the 
Cotton-Mouton effect measurement in helium.— We will 
however detail the main steps, since a slightly different 
method was used in the present case. 

To extract the ellipticity 4'cM(t) from Eq. (fT0)l . we cal¬ 
culate the following Y{t) function 


Y{t) 


A,t(t) 


-h 


DC 


2|r| 


7'f'CM(i) -I- 


^cm(0 

2|r| 


|e|eF(t) 

^ 2|r| 


e|(t) 

2|r| ’ 

( 20 ) 


A. Magnetic field 

The transverse magnetic held B± is generated by an 
X-Coil, specially designed by the High Magnetic Field 


where 7 stands for the sign of T. Jdc is the static sig¬ 
nal measured just before the application of the magnetic 
held. The absolute value of the static ellipticity |r| is 
also measured before each pulse. 
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Two parameters are adjustable in the experiment: the 
sign 7 of the static ellipticity T and the direction of the 
transverse magnetic field. We acquire signals for both 
signs of r and both directions of Bj^: parallel to x is 
denoted as > 0 and antiparallel is denoted as < 0. This 
gives four data series: (F > 0, B± > 0), (F > 0, Bj_ < 0), 
(F < 0, B± < 0) and (F < 0, B± > 0). 

For each series, signals calculated with Eq. (l20ll are av¬ 
eraged and denoted as !>>, T><, F<< and To- The 
first subscript corresponds to F > 0 or < 0 while the sec¬ 
ond one corresponds to parallel or antiparallel to x.. 
This average function can be written in a more general 
form than the one of Eq. (PIB . It is the sum of different 
effects with different symmetries, denoted as s 


r» = +* + !( 

= -» + !( j^)»++ 



( 21 ) 


The first subscript in s corresponds to the symmetry with 
respect to the sign of F and the second one to the symme¬ 
try with respect to the direction of Bj^. The subscript 
-|- indicates an even parity while the subscript — indi¬ 
cates odd parity. The ratio < 1/F > is the average of 
1/|F| measured during corresponding series. The terms 
and 0p are included in s++, 7|e|0F are included 

in s _, and s-|_corresponds to a spurious signal with 

an odd parity towards the direction of Bj^ and an even 
parity with respect to the sign of F. The ellipticity 7d>cM 
corresponds to s _ 

From this set of four equations with four unknown 

quantities (4 'cm, S-i-i-, s _and S-|_), we extract 4'cM(t), 

which is fitted by aB\ f. The cavity filtering should 
again be taken into account, as indicated by the subscript 
f The Cotton-Mouton constant fccM finally depends 
on the measured experimental parameters as follows: 


TABLE III. Parameters that have to be measured to infer the 
value of the Cotton-Mouton constant fccM and their respective 
relative A- and B-type uncertainties at la. Typical values are 
given at P = 8 X 10“® atm. 


Parameter 

Typical 

value 

Relative 

A-type 

uncertainty 

Relative 

B-type 

uncertainty 

r [ms] 

1.14 

2.0 X 10"^ 


a X 10® [T-2] 

2.82 

2.8 X lO”'^ 

2.2 X 10"^ 

[MHz] 

65.996 


3 X 10"'^ 

Lb [m] 

0.137 


2.2 X 10"^ 

A [nm] 

1064.0 


< 5 X 10"'‘ 

sin 20p 

1.0000 


9 X lO""* 


kcM 2.31 

2.0 X 10"^ 

3.1 X 10"^ 

^IQie [rp-2j 




TABLE IV. Parameters and their respective relative A- and 
B-type uncertainties at la that have to be measured to infer 
the value of the normalized Cotton-Mouton constant 


Parameter 

Typical 

value 

Relative 

A-type 

uncertainty 

Relative 

B-type 

uncertainty 

kcM 

xiQie [7-2j 

2.31 

2.0 X 10“^ 

3.1 X 10"^ 

P X 10® 

[atm] 

5 


2 X 10"® 

linear fit 
xl0®'‘ 

[T-^atm"®] 

2.41 

1.5 X 10"® 


^CM 

2.41 

1.5 X 10"® 

3.1 X 10"^ 


[T-^atm-i] 


pressure are fitted by a linear equation, and we obtain for 
the value of the Cotton-Mouton constant at P = 1 atm 


fccM(r,P)= , afsr A^^T- (22) 
47rrA^°^ Lb sm20p 

C. Measurement and error budget 

The A- and B-type uncertainties associated to the mea¬ 
surement of kcM are detailed in Tab. lIIIl and are given at 
1(7. The B-type uncertainties have been evaluated previ¬ 
ously and detailed in Ref. [i3- They essentially come from 
the length of the magnetic field Lb and the fit constant 
a. 

We have measured the Cotton-Mouton constant in 
xenon at T = (293 ± 1) K and for nine pressures from 
3 X 10“^ to 8 X 10“^ atm. The data as a function of the 


= (2.41 ± 0.37) X 10”^^ T-^atm-^ (23) 

The uncertainty given at Icr is detailed in Tab.|TVl The 
dominant uncertainty comes from the linear fit of the 
Cotton-Mouton constant versus pressure (A-type). The 
value of fcgjy; normalized at 273.15 K is calculated with a 
scale law on the gas density 

= (2.59 ± 0.40) X 10”^^ T-^atm-^ (24) 


V. OUR CALCULATIONS 

The Verdet constant and the Cotton-Mouton birefrin¬ 
gence were computed within Coupled Cluster response 
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theory^ at the CCSDl-^ and levels of ap¬ 

proximation. Specifically, the Verdet constant was ob¬ 
tained from the following frequency-dependent quadratic 
response function^i^i^— 


For further details on how the above Cauchy moments 
and dispersion coefficients of the given quadratic response 
function are computed within coupled cluster response 
theory, the reader should refer to Refs. 123 - 1^ . 


= (25) 

with C = „ — = 0.912742 x 10“^ in atomic units, 

N the number density [N = for ideal gases), e the 
elementary charge, me the electron mass, cq the speed of 
light in vacuo, the frequency of the probing light, 

and [ix,y and are Cartesian components of the elec¬ 
tric dipole, and angular momentum operators, respec¬ 
tively. The hypermagnetizability anisotropy A 77 enter¬ 
ing the Cotton-Mouton birefringence in Eq. (|4]) (the only 
term contributing for atoms) is given by the combination 
of a quadratic and a cubic response functions^ 


= A77P-h 


^ {{l-^x'i l-txj ^xx))ui,0 

(26) 


with Qxx the Cartesian component of the traceless 
quadrupole operator. At the CCS level, calculations were 
performed at three different wavelengths, namely 1064, 
632.8 and 514.5 nm. At the CCSD level, we computed the 
dispersion coefficients, as done in our previous studypi 
i.e., for the Verdet constant 

V{2n) =2nS(-2n-2); (27) 

00 

V(a;) =C'^a;2"V(2n); (28) 

n—1 

whereas for the Cotton-Mouton constant 
AT](2n) = —-[(2n -|- l)(2n -|- 2)S{—2n — 4) -|- B(2n)]] 

(29) 

00 

Ar]{uj) = '^uj'^'^Ar]{2n). (30) 

n =0 


Relativistic effects were approximately accounted 
for by employing relativistic effective core poten¬ 
tials (ECPs)^! and specifically pseudo-potentials (PP). 
“Small core” effective pseudo-potentials were used to de¬ 
scribe the 28 inner electrons (that is, the [Ar]3d^° core), 
whereas the remaining 26 valence electrons were corre¬ 
lated as in standard non-relativistic calculations. The 
basis sets used were constructed starting from the singly 
augmented aug_cc_pvxz_pp (x=t,q) sets of Peterson et 
al^ Since single augmentation is usually not sufficient to 
ensure converged results, at least for the Cotton-Mouton 
birefringence. Additional sets of diffuse functions were 
added by applying an even-tempered generation formula 
commonly used for this purpose to the orbital func¬ 
tions describing the valence electrons, while retaining the 
pseudo-potential of the original set. The resulting sets 
are labeled d-aug and t-aug, for double and triple aug¬ 
mentation, respectively. 


Where pseudo-potentials parametrically account for 
relativistic effects on the innermost orbitals, other rel¬ 
ativistic effects (e.g. higher-order and picture change 
effects, spin-orbit coupling) could play a significant 
rolej ^^i^° When dealing with valence properties like elec¬ 
tric hyperpolarizabilities, the higher-order relativistic ef¬ 
fects and picture change effects (for the dipole operator 
and also the electron-electron interaction) are expected 
to be not so important. Also, spin-orbit coupling should 
be quite weak. Both the Faraday and Cotton Mouton 
birefringences, however, involve the magnetic dipole op¬ 
erator. In general relativistic effects on magnetic prop¬ 
erties can be more significant and more difficult in terms 
of picture change (the operators look different in rela¬ 
tivistic and non-relativistic theory and this may require 
a correction of the property operator that one uses as a 
perturbation) 


Above, S{k) is the Cauchy moment 

Sik)=Y^ 2a;^’o^(0 \ fi, \ m){m \ \ 0) (31) 

m /0 

with huJmo indicating the excitation energy from the 
ground state 0 to the excited state(s) m, and B{2n) is 
the dispersion coefficient introduced when expanding, for 
frequencies below the lowest excitation energy, the elec¬ 
tric dipole-electric dipole-electric quadrupole quadratic 
response function Rx.x,xx(-w; w, 0 ) = {{fix', ^J^x,Qxx))uJ,o 
in a convergent power series in the circular frequency ui 


Nonetheless, also given that the most stringent require¬ 
ment in terms of basis set convergence is the inclusion of 
diffuse functions as in the case of the electric hyperpo- 
larisability, it is reasonable to assume that both proper¬ 
ties are essentially valence properties, for which picture 
change effects are typically small, and we reckon there¬ 
fore that the use of (PP)ECPs can be considered accurate 
enough. 

The results obtained in the x=q basis sets are summa¬ 
rized in Tab. Eland Tab. ED for CCSD and CC3, respec¬ 
tively. 


Bx,x,xxi-uj;uj,0) = y^cu^"R(2n) 

n =0 


(32) 


All calculations were performed with the Dalton 
code<^ 



7 


TABLE V. Dispersion coefficients of the Verdet and Cotton-Mouton response functions at the CCSD level of theory (atomic 
units). 


n 

B{2n) 

S'(-2n-4) 

S{-2n - 2) 

V{2n) 

AT](2n) 

aug cc pvqz pp 

0 

-654.89471 

126.50595 



100.47070 

1 

-8903.3825 

763.59899 

126.50595 

253.01190 

-64.951345 

2 

-92298.251 

5369.0486 

763.59899 

3054.3960 

-17193.302 

3 

-860869.9 

41692.560 

5369.0486 

32214.292 

-368478.36 

4 



41692.560 

333540.48 


d-aug cc pvqz pp 

0 

-739.15630 

126.97174 



121.30323 

1 

-9822.9127 

774.87190 

126.97174 

253.94348 

131.11247 

2 

-106369.19 

5553.2321 

774.87190 

3099.4876 

-15056.943 

3 

-1074975.1 

44095.369 

5553.2321 

33319.393 

-348591.39 

4 



44095.369 

352762.95 


t-aug cc pvqz pp 

0 

-748.34187 

126.91927 



123.62583 

1 

-9940.7218 

774.47234 

126.91927 

253.83854 

161.76343 

2 

-107513.63 

5551.3771 

774.47234 

3097.8894 

-14756.921 

3 

-1084127.2 

44088.280 

5551.3771 

33308.263 

-346204.12 

4 



44088.280 

352706.24 



TABLE VI. CCS values of the 

response function components (in atomic units) involved in 

the Verdet and Cotton-Mouton 

birefringences. The Verdet constant V”(ta) is 

given in atm ^ rad. T’ 

^m ^ and the Cotton-Mouton constant 

7 N • • 

ftcM IS in 

T ^.atm 

1 at 273.15 K. 







A[nm] 

l^y-) Lz))uj,Q 

V”(ca) X 10^ 

l-^xi ^xx)^tjj,0 

((/^X) Lzi Lz')')uj,U. 

;,0 

At; k 

N 

CM ^ 

aug cc pvqz pp 

1064 

11.1587 

3.505 

-668.242 

272.564 


98.9195 

2.239 

632.8 

19.5823 

10.35 

-700.706 

308.069 


98.1593 

2.222 

514.5 

24.9438 

16.22 

-728.260 

339.617 


97.1607 

2.200 

d-aug cc pvqz pp 

1064 

11.2155 

3.522 

-755.936 

274.099 


120.459 

2.727 

632.8 

19.6927 

10.40 

-791.994 

310.285 


120.427 

2.726 

514.5 

25.0963 

16.32 

-822.705 

342.514 


120.048 

2.718 

t-aug cc pvqz pp 

1064 

11.2127 

3.521 

-765.680 

274.031 


122.912 

2.782 

632.8 

19.6878 

10.40 

-802.186 

310.210 


122.994 

2.784 

514.5 

25.0901 

16.31 

-833.274 

342.435 


122.710 

2.778 


VI. RESULTS AND DISCUSSION 
A. Faraday effect 
1. Experiments 

We can compare our value of the normalized Verdet 
constant to other published values. The most extensive 
experimental compilation of Verdet constants has been 
reported by Ingersoll and Liebenberg in 1956, for sev¬ 
eral gases including xenon^^ for wavelengths ranging from 
363.5 to 987.5 nm, with a total uncertainty of about 1 %. 
These values are plotted in Fig.O 

No datum has ever been reported for A = 1064 nm. 
Nevertheless, we can extrapolate its value from the points 


of Fig.[2l by fitting the data with a function of form 
V = Aj}? + BjX^ (solid curve in Fig.[2))i^‘^ A sup¬ 
plementary systematic uncertainty should also be added, 
since the authors measured the ratio between Faraday 
effects in xenon and in distilled water, and rescaled their 
measurements with accepted values for water . Thus 
it does not correspond to absolute measurements of the 
Faraday effect, contrary to ours. 

At A = 1064nm and T = 273.15K we obtain = 
(3.46 ± 0.04) X 10“^ atm“^rad.T“^m“^. The Icr uncer¬ 
tainty includes the one given by the fit. This value is 
compatible with our experimental value (Eq. (11811 ). rep¬ 
resented as the open circle in Fig. [5] and as the straight 
and dashed lines in Fig. [31 























FIG. 2. A: Experimental values of xenon normalized 
Verdet constant at T = 273.15 K reported by Ingersoll and 
Liebenberg^^ for wavelength from 363 nm to 987.5 nm. These 
values are htted by the law Aj}? + B (solid line), o: Our 
experimental value at T — 273.15 K. 


E 



i Savukov 



3.8x10“^ - 




3.7 - 




3.6 - 


Ikalainen et al 


3.5 - 

Ekstrom et al 

• NR-CCSD 
»X2C 

• CCS 

• CCSD 

3.4 - 

TDDHF 

DHF 

Thi.s work 


TDHF* 

• NR 



TABLE VII. Experimental and theoretical values of the nor¬ 
malized Verdet constant at T = 273.15 K, A = 1064 nm, with 
uncertainties at Icr. 


Ref. 

V" X 10^ 

Remarks 


(atm~^rad. 

T-V-i) 



Experiment 


Ingersoll et 3.46 ± 0.04 Interpolated with 

AjX^ + 2B/X*. 
Scaled to water. 

This work _ 3.56 ± 0.10 _ 

Theory 


Savuko\^ 

3.86 ±0.01 

Interpolated in this work 
with A/A^ ±R/A'‘. 

Ekstrom et ai^ 

3.35 

TDHF 

Ekstrom et al^ 

3.46 

TDDHF 

Ikalainen et al^ 

3.34 

NR 

Ikalainen et al^ 

3.48 

X2C 

Ikalainen et al^ 

3.46 

DHF 

Ikalainen et al^ 

3.52 

NR-CCSD 

This work 

3.49 

C C S D /t- aug_cc_pvqz _pp 

This work 

3.52 

CC3/t-aug_cc_pvqz_pp 


FIG. 3. Normalized Verdet constant of xenon at T = 273.15 K 
at A = 1064 m. Solid line: our experimental mean value. 
Dashed lines: our experimental value with Icr uncertainty. 
Points : theoretical predictions (both ours and from the lit¬ 
erature). See text and Tab. lVIll for the references. 


2. Theory 

We can also compare our experimental value with the¬ 
oretical predictions (both ours and from the literature), 
plotted in Fig.|3] and summarized in Tab. lVIll at 1 atm, 
273.15 K and with the gas number density of an ideal 
gas. To convert from theoretical results given in atomic 
units into the units used experimentally, we exploited the 
relation: 

V (atm-^rad.T-V-i) = V (a.u.) x 8.039617 x 10^ 

(33) 

Our experimental value is compatible within Icr with 
both our “best” coupled cluster results (t-aug_cc_pvqz_pp 
basis) and the theoretical prediction of Ekstrom et al^ 
within 2(7 with the estimate of Ikalainen et al^ and 
within 3 ct with that of Savukovi^ 

The uncertainty of a few percent obtained on our ex¬ 
perimental value allows to comment on the agreement 
with theoretical predictions as a function of the the¬ 
oretical approximation or model. Savukov^^ has used 
a relativistic particle-hole configuration interaction (Cl) 


method. He does not give a value at 1064 nm, but the 
latter can be interpolated, as done with the previous ex¬ 
perimental data of Ingersoll and Liebenberg)^ obtain¬ 
ing the value of Tab. Em with an uncertainty given by 
the fit. The agreement between theory and experiment 
is only within 3cr, even if relativistic effects are taken 
into account. Ekstrom et al^ have used the nonrela- 
tivistic time-dependent Hartree-Eock (TDHF in Fig.|3]) 
and the relativistic time-dependent Dirac-Hartree-Fock 
(TDDHF in Fig. [3]). There is clearly a better agreement 
(better than Icr), between their calculations and our ex¬ 
perimental value when relativistic effects are taken into 
account. Finally, Ikalainen et al^ have used the non- 
relativistic Hartree-Fock method (NR in Fig. [3]), the ex¬ 
act two-component method (X2C in Fig.|3]), and the fully 
relativistic four-component method (DHF in Fig. [3]). The 
same authors also report (in the supporting information 
file) a non relativistic CCSD result (NR-CCSD in Fig.[3|). 
While their uncorrelated results confirm that relativistic 
effects should be taken into account to improve agree¬ 
ment with experiment, their non-relativistic CCSD re¬ 
sult highlights how the inclusion of correlation effects is 
equally important. Also worth noticing is the rather poor 
performance of the BLYP and B3LYP functionals, which 
overestimate the value of the Verdet constant in both 
non-relativistic and relativistic calculations. This also 
applies for the BHandHLYP functional in the relativistic 
calculations, whereas the non-relativistic BHandHLYP 
value is still within Icr of our experimental result (See 
Table S5 of the Supporting Information file of Ref. l34) . 
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TABLE VIII. Experimental (uncertainties of la) and the¬ 
oretical values of the Cotton-Mouton constant of xenon at 
T = 273.15 K. 


Ref. 


A (nm) 

k^CM X 10"" 
(T-^atm-i) 

Experiment 

Carusotto et al^ 


514.5 

(2.29 ±0.10) 

Hiittnep^ 


632.8 

(2.41 ±0.12) 

Bregant et 


1064 

(3.02 ±0.27) 

This work 


1064 

(2.59 ± 0.40) 

Theory 

Bishop et ai^ 


OO 

2.665 

This work, 


514.5 

2.803 

CCSD/t-aug_cc_pvqz_pp 
This work, 


632.8 

2.808 

CCSD/t-aug_cc_pvqz_pp 
This work, 


1064 

2.804 

CCSD/t-aug_cc_pvqz_pp 
This 

work, 

514.5 

2.778 

CC3/t-aug_cc_pvqz_pp 

This 

work, 

632.8 

2.784 

CC3/t-aug_cc_pvqz_pp 

This 

work, 

1064 

2.782 

CC3/t-aug cc pvqz pp 





B. Cotton Mouton Effect 
1. Experiments 

Only a few measurements of the Cotton-Mouton effect 
in xenon have been discussed in the literature. There is 
one at A = 514.5 nm by Carusotto et al^ one at A = 
632.8 nm by Hirttner (reported as a private communica¬ 
tion by Bishop et al)^ and finally one at A = 1064 nm 
by Bregant et Our experimental value, referring 

to A = 1064 nm is compatible within la with the datum 
of Refs. The set of results is shown in Tab. lVIIll 

and plotted as a function of the wavelength in Fig.lH 

Our measurement has an uncertainty of about 15%. 
This value, which is larger than that of the other re¬ 
ported values, especially those given for wavelengths of 
514.5 nm and 632.8 nm, was established via a complete 
error budget. Note that no information is available on 
the setup, the number of pressures, the error budget and 
the evaluation of the uncertainty for the value reported 
at A = 632.8 nm by Bishop et as a private commu¬ 
nication of Hiittner. The value reported at A = 514.5 nm 
by Carusotto et al^ was measured only at 1 atm, and 
by comparing the observed magnetic birefringence with 
that of nitrogen under the same experimental conditions, 
therefore taking as a reference, free of uncertainty, the 
Cotton-Mouton constant of nitrogen. It is safe to say 
therefore that the uncertainty associated to their datum 
might be underestimated. Finally, the value reported by 



X (nm) 

FIG. 4. Reported values of Cotton-Mouton constant of xenon 
for A ranging from 514.5 nm to 1064 nm and with Icr un¬ 
certainty. Experimental values: black triangle: Carusotto 
et al^ open triangle: Hiittner (private communication by 
Bishop et al.)^ black diamond: Bregant et open di¬ 

amond: this work. Theoretical predictions: dashed line: SCF 
method for A = oo by Bishop,— open circle: this work, CCSD, 
black circle: this work, CCS 


Bregant et at A = 1064 nm corresponds to the 

weighted average between measurements at two differ¬ 
ent pressures (9 pressures for our measurement) and the 
uncertainty is similar to ours. 


2. Theory 

The Cotton-Mouton constant fccM is linked to Ary by 
the relationship^ 

, , -irT.-2x 6.18381 x 10-14^ , , 

fccM (atm It = - - - Ar] (a.u.). (34) 

Only one theoretical prediction has been published so- 
far for the Cotton-Mouton effect in xenon— The cal¬ 
culation of Bishop and Cybulski was performed at the 
self-consistent-field (SCF) level of approximation, and it 
yielded static hypermagnetizability anisotropy Arj. As 
stated by the authors, relativistic effects were not taken 
into account, even though the authors expected them to 
play a substantial role. Our experimental value agrees 
with that theoretical prediction within la. 

Our computed coupled cluster results, both CCSD and 
CC3, in the largest (t-aug_cc_pvqz_pp) basis sets for the 
three wavelengths at which experimental results are given 
in Tab. IVIIII Both the CCSD and CC3 values at 1064 
nm are well within Icr of our experimental measurement, 
and just outside Icr of the result by Bregant et al.—’^ 
At 632.8 nm the agreement of our CC3 value with the 
experimental result of Hiittner— is just outside 3cr. At 
514.5 nm our computed values fall well outside 3cr of the 
estimate of Carusotto et al^ This apparently confirms 
that the error associated to this measured value might 
be underestimated. 
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VII. CONCLUSION 

We have carried out a thorough analysis of the Faraday 
(circular) and Cotton Mouton (linear) birefringences of 
xenon, at a wavelength of 1064 nm. The study involves 
both an experimental segment, exploiting the capabilities 
of a state-of-the-art optical setup, and a computational 
element, where sophisticated wavefunction structure and 
optical response models (and with an estimate of the ef¬ 
fect of relativity) were employed. 

Our experimental estimate for the normalized Verdet 
constant of xenon at a temperature of 273.15 K and 
A=1064 nm, 1/^ = (3.56±0.10) xlO"^ atm"^ rad 
m“^, is very well reproduced by our theoretical approach, 
which yields a value = 3.52 xl0“^ atm“^ rad 
m“^ using the CC3 approximation) within la of the mea¬ 
sured datum. 

With respect to the Cotton Mouton effect, at 
T=273.15 K and A=1064 nm, experiment yields a nor¬ 
malized constant = (2.59zb0.40) xl0~^^ atm~^ T“^, 
whereas we compute (again with our most sophisticated 
model, CC3) a value of = 2.78 xl0“^'^ atm“^ T“^, 
therefore within la of experiment. 
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